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SUMMARY 

A method is described for determining the "displacement surface" 
of a known three-dimensional compressible boundary- layer flow in terms 
of the mass-flow defects associated with the profiles of the two 
velocity components parallel to the surface . 

The height of the displacement surface above the body surface for 
flow about a yawed infinite cylinder is shown to be equal to the height 
characterizing the mass-flow defect of the chordwise velocity profile. 
The displacement surface height is shown to differ, in general, from 
that associated with the resultant mass-flow defect, even at stagnation 
points of the secondary flow. Numerical values are found for the known 
three-dimensional houndary-layer flow about a cone at a small angle of 
attack to a supersonic stream. 


INTRODUCTION 

The boundary layer established in the flow of a slightly viscous 
fluid about a body is normally considered an Isolated region -therein 
the effects of viscosity predominate, and. outside of which the motion 
of the fluid is governed by the laws of nonviscous motion. For large 
Reynolds number, the boundary layer is assumed to be so thin that the 
nonviscous portion of the flow occurs as though there were no boundary 
layer. This assumption is strictly correct in the limit of infinite 
Reynolds nufiiber. For large hut finite Reynolds numbers, the growth of 
the boundary layer causes the stream to be deflected away from the body 
surface. 

This displacement effect of the boundary layer on the nonviscous 
flow may properly he determined from the behavior of the boundary layer 
itself, as established either by eaperiment or by solution of the 
Prandtl boundary-layer equations for laminar flow. 
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It does not follow, however, that this revised outer flow may he 
used to improve the solution for the boundary layer, still using the 
Prandtl equations. A new set of equations must he used for this pur- 
pose t akin g into account the variation of pressure across the boundary 
layer, which is neglected in the Prandtl equations. (See Alien* s iter- 
ative solution for incompressible flat-plate flow, reference 1. ) 

The customary definition of displacement thickness (reference 2) 
is applicable to two-dimensional flow and is expressed in the following 
equation: 



pu dy = p 1 u 1 (h-&*) 


( 1 ) 


where h is some .location well outside' the boundary layer (fig. l), 
at which pu = Pi u i, and beyond which the flow may be considered non- 

viscous. (A list of symbols is provided in the appendix. ) Under the 
assumption of an extremely thin boundary layer, h is so Bmall that 
PjUj may be taken as the evaluation at the body surface (y=0) of the 

nonviscous flow obtained by neglecting the presence of the boundary 
layer. Equation (l) equates the actual mass flow near the surface with 
the mass flow which would be associated with a nonviscous flow that 
terminates at 5* rather than extending to the wall. Thus, the non- 
viscous portion of the flow behaves as if it occurred in the presence 
of a solid boundary given by the displacement thickness 5*(x) . Equa- 
tion (l) may be solved for 5*: 



Ordinarily, theoretical boundary-layer solutions for pu/ p^Up asymp- 
totically approach 1 for large values of Reynolds number based on y. 
Therefore, displacement thickness is often defined as follows: 
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For three-dimensional h oundary- layer flows, two lengths character- 
izing mass-flow defect may he defined in terms of the profiles of the 
two velocity components tangential to the surface (fig. 2), 



DO 

o) and it is not clear -which, if either, defines a displacement surface 

that properly describes the extent to which the nonviscous flow is 
deflected by the boundary layer. 

Of course, it is expected that there does exist a displacement 
surface for such flows. The analysis that follows shows that such a 
surface may be described, using a defining equation more fundamental 
than equation (l). 

The velocity and density profiles are assumed to be known for the 
three-dimensional flows under consideration. Qases for which this is 
true include: the laminar boundary layer on yawed infinite cylinders, 

treated by Prandtl (reference 3), R. T. Jones (reference 4), and by 
Sears (reference 5); the laminar boundary layer on a cone at a mnnii 
angle of attack to a supersonic stream (reference 6). The displacement 
effect of these flows will be treated specifically. 

This report may be considered as an addendum to references 6 aurl 7. 
This investigation was conducted at. the HACA Lewis laboratory. 


THEORT 

Defining Relation for Displacement Surface 

The boundary- layer solution (assumed known) yields a certain dis- 
tribution of velocity v^ ^ (x, z ) normal to the body surface at the 

outer edge h(x, z) of the boundary layer, where p, u, and w may be 
taken essentially equal to p 1 , u^, and w 1 (see fig. 2). Under the 

assumption that the nonviscous flow is altered only slightly by dis- 
placement, the most direct way to compute this effect would be to sup- 
pose that the nonviscous equations hold for y5jh(x, z) and to impose 
the following boundary condition on the normal velocity v 0 p in the 
outer flow: 


v of * T b i( x > z ) at y = Mx,z) 


( 3 ) 
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Hcweve^ since the boundary conditions usua ll y encountered in non- 
viscous flow specify an impermeable surface, it is convenient to recast 
the boundary condition given in equation (3) in answer to the question: 
What impermeable surface y = A(x, z) would deflect a nonviscous fluid 
in such a way as to produce a normal velocity satisfying condition (3)? 
This fictitious surface may be called the displacement surface. 

Because A(x, z) is imagined to be a stream surface, it is neces- 
sary to specify that at y=A the resultant velocity vector (u,v,w) 
be tangent to the surface y=A. Thus, at y=A the ratio of v to the 
magnitude of _q, which is defined as the vector (u,w), must be set 
equal to the slope of the surface y=A, measured in the direction of 
the vector jj, or, equivalently, equal, to the component in the direction 
of ^ of the vector grad A. In vector notation, therefore, the nor- 
mal velocity v which would be produced in a nonviscous fluid at an 
impermeable surface y=A is . _£ • grad A. The vector q may be 
obtained by evaluating the velocity vector of the unrevised nonviscous 
flow at the body surface y=0, under the related assumptions that the 
velocity vector varies only slightly over distances of the order of 
the actual boundary-layer thickness, and that the revision required to 
take account of displacement is slight. Thus, at y=A (see fig. 3), 

v Qf - * grad A 

The increment in v 0 f between A and h is approximately (h-A) 5v of /Sy, 

again assuming a thin boundary layer and using only the first term in a 
Taylor series. 

# 

To the order of approximation contemplated in this analysis, 
dv of /dy may be obtained from the unrevised nonviscous flow evaluated 

at y=0. Thus, the fictitious impermeable surface A would produce, 
at y=h, a normal velocity (see fig. 3) 


of 


= % 


• grad A + (h-A) 



Introducing this result into the boundary condition given in equa- 
tion (3) yields the defining relation for A(x, z): 


(l, bf>y=h 



+ _3]_ ' grad A 


( 4 ) 


The boundary- layer solution yields v^j. 
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Expression of A in Terms of Mass-Flow Defects 


The displacement surface A may be related to the mass-flow 
defects (equations (2)) which characterize the boundary layer as 
follows: in a Cartesian coordinate system (fig. 2) , the equation of 

continuity for both the boundary- layer and nonviscous flow is 


Spv Spu Spw 

W = ^ ~ ST 


(5) 


Under the Prandtl boundary-layer assumptions, the Cartesian equations 
of motion, and hence equation (5), may be applied in an orthogonal 
curvilinear coordinate system in which the surface of the body is 
given by y=0, provided the radius of curvature of the body is every- 
where large as compared with the boundary-layer thickness. Integrating 
equation (5) across the boundary layer yields 





dy 



or, inasmuch as h Is outside the boundary layer where pu = Pp u l and. 
pw = PqW-l, 





S 


(Pl w l 5 z) 


( 6 ) 


where S x and 6 Z are the mass-flow defects defined in equations (2). 

For the nonviscous flow, v D f =0 at the body surface (y=0) , and equa- 
tion (5) becomes 



f dPl u l 


§z~ ) 


(?) 
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Introducing equations (6) and (7) into equation (4) yields 
[Pl-1(A-5 X )] + ^ [p 1 v 1 (A-B z jj = 0 


In a Cartesian coordinate system, the displacement surface is related 
to the mass-flow defects hy equation (8a). 'When cases arise for which 
other coordinate systems must he used, the following generalization in 
■vector form of equation (8a) may he used: 


div 



PiSi A - / (piSi-p.2) & 


0 


(at) 


where y is the distance normal to the hody surface and where the div- 
ergence operator involves differentiation only with respect to the two 
coordinates parallel to the hody surface. 


Examples 


Plane flow . - Equation (8a) may he integrated to yield , for plane 
flow (d/dz = 0) : 


A = 


8 X + 


K 


Pl u l 


( 9 ) 


where K is a constant of Integration. 

The appearance of this constant means, in general, that the 
revised houndary condition on the nonviscous flow near the wall may he 
applied at any surface in the houndary-layer region; for example, 
along the wall Itself. Of course, if there is a stagnation point on 
the hody where u^ vanishes, then K must he taken equal to zero and 

the revised houndary condition must he imposed at the location A = 5^, 
at least near the stagnation point. 

Stagnation point of secondary flow. - In the vertical plane of sym- 
metry of the flow about a hody of revolution at an angle of attack, the 
circumferential velocity component vanishes in the boundary layer as 
well as in the nonviscous outer region. Such stagnation of the second- 
ary motion would occur in a variety of cases, in particular, wherever 
an essentially three-dimensional flow has a plane of symmetry. "When 
Wj_ is designated as the component of secondary flow (circumferential 

velocity for a hody of revolution) at the outer edge of the boundary 
layer, and equation (8a) is written in the form 


2536 
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^ [PXUX^)] + "1 ^ [pl^Sj] + Pl(A-S z ) |S = 0 

it is clear that A cannot equal V S X , in general, even if both w^_ and 
the boundary- layer profile of w vanish, unless cSwjJdz also vanishes. 

Flow about yawed infinite cylinder . - In this case, there is a 
spanwise boundary- layer velocity profile and an associated spanwise 
mass-flow defect. (See references 3 to 5. ) If x is taken in bhe 
chordwise direction (fig. 4), the entire flow depends only on x. Thus, 
derivatives with respect to the spanwise coordinate z vanish) and, 
hence, from equation (8a), the p].ane -flow result (equation (9)) applies. 
Accordingly, the spanwise mass-flow defect represented by 6 Z does not 

enter into the determination of the displacement surface. 



Because the outer flow is conical, p^., u^, and w^ are functions only 
of cp) and equation (ll) may be written 


sin 0 [x(A-S x )J + JL ^^(A-Stpjj- = 0 (13) 

In the case of conical outer flow, the associated boundary- layer 
profiles show similarity of the Blasius type in meridional planes (see 
reference 7 or 8). Thus, in a meridional plane. A, 5 x , and 6 are 

proportional to tfx. Incorporation of this information into equa- 
tion (l3) gives 

| Pqu-L sin 0 (A-5 X ) + ^ ^pjwjA-b^J = 0 


(14) 
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For a cone at; small angle of attack a, up is nearly equal to u, 
the velocity on the cone surface at zero angle of attack. The quanti- 
ties p^, A, and 6^, vary only slightly with angle of attack, 

whereas 


Wq_ = ecu Ag sin <p 


(15) 


where A£ depends only on the cone vertex angle and the flow Mach num- 
ber and is defined in reference 6 as follows: 


c u u sin © 

The quantities z, x, and u are in the notation of reference 9, 
wherein they are tabulated. 

To the first order in a, the substitution of equation (15) into 
equation (14) yields the result 


A = &x + | A 2( V S x) COG 'P ( 16 ) 

The analysis qf reference 6 yields the values of 5 X and 8<p. 

Clearly, A differs from 8 X in the plane of symmetry <p = 0, it, 
where the circumferential velocity w vanishes. It might, however, be 
noted that in the plane T = -g, — , A = 8 X . Figure 6 shows the pro- 
portional rate of increase of displacement thickness A with angle of 
attack in the plane of symmetry <p = it at zero angle of attack. The 
corresponding rate of increase of the mass-flow defect 8^. is shown 

for comparison. These curves are obtained from equation (16) and the 
results of reference 6. Hie sine of the semivertex angle is introduced 
as a factor primarily to permit presentation of the case 0 = 0 as the 
limit of an indeterminate form. For a stream Mach number of 2, the 
change in displacement thickness appears to be of the order of 50 per- 
cent larger than the change in the meridional, mass-flow defect. 


COXICIUDOTG- REMARKS 

The foregoing analysis deals only with the displacement effect of 
a known boundary layer on the nonviscous outer flow, and hence can be 
applied only if the boundary-layer behavior has been determined either 
theoretically or experimentally. The latter approach might possibly 
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find application in the correction of nozzle contours for boundary- layer 
development. In this connection, perhaps, it should he noted that the 
analysis is not restricted to 1 ami nn.r flows. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 6, 1952 
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APPENDIX - SYMBOLS 

The following symbols are used in this report: 

Ag function of cone angle and Mach number (equation (15)) 

h height above wall at which p, u, w = p^, u-]_, w-j_, and beyond which 

nonviscous equations apply (equation (3)) 

Mq stream. Mach number 

^ velocity vector composed of components parallel to body surface 
u,w 

u velocity component in x-direction 

u meridional velocity component at surface of cone at zero angle of 
attack 

v velocity component in y-d l rection 

vq,z boundary- layer solution for velocity normal to surface , evaluated 
at outer edge of boundary layer 

v 0 f nonviscous solution for velocity normal to surface, evaluated near 
the surface 

w velocity component in z-direction 

X \ 

/ coordinates in body surface 

z J 

y coordinate normal to surface 

a angle of attack 

A height above body of displacement surface 

6* displacement thickness in plane flow 

5 length characterizing mass-flow defect of u-profile (equa- 
tions (2), (12)) 





2536 


MCA TN 2722 


11 


5 Z length characterizing mass-flow defect of w-profile 

(equation (2)) 

5<p length characterizing mass-flow defect of w-profile 
(equation (12)) 

0 semivertex angle of cone 

p density 

<p angular .coordinate (fig. 5) 

Subscript 1 denotes evaluation of nonviscous flow at body surface, 
taken equivalent to conditions at outer edge of a boundary layer of 
infinitesimal thickness. 
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Figure 1. - Plane boundary- layer flow. 









Figure 4. 


Coordinate system for yawed in_ 


> 





Figure 5. - Coordinate system for cone at angle of attack 





